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Abstract -In this Letter we report how thermodynamic properties of a giant frustrated magnetic 
Keplerate molecule of A'^ = 30 spins s = 1/2 can be evaluated with the help of the highly accu- 
rate finite-temperature Lanczos method. The comparison to experimental data shows excellent 
agreement. Since this molecule is structurally related to the archetypical kagome lattice anti- 
ferromagnet we expect new detailed insight into properties of this important class of frustrated 
materials. 



Introduction. — Nanometer sized polyoxometalate 
molecules constitute a fascinating class of molecular mate- 
I rials (l}|7j. The series of Keplerate molecules {Mo72Fe3o}, 
{Mo72Cr3o}, {M072V30}, {W72V30} is from a magnetism 
^ point of view of special interest since in these bodies para- 
Q magnetic ions occupy the vertices of a nearly perfect icosi- 
I— dodecahedron - one of the Archemidean solids. These 
^ bodies resemble some of the most interesting magnetically 
>^ frustrated spin lattices such as the kagome lattice antifer- 
romagnet [8]. Valuable insight about the physics of such 
00 lattices can be gained by studying the finite-size bodies. 
O^But although icosidodecahedra consist of only = 30 
'^spin sites, the dimension of the Hilbert space for s — 1/2 
CN reaches a stunning 1,073,741,824. 

" } Figure [T] shows the structure of the icosidodecahedron: 
^^spin sites are displayed by bullets, edges are given as 
• -straight lines - in the later used Heisenberg model they 
p'j represent the interaction pathways. If such interactions 
are of antiferromagnetic nature, i.e. favour antiparallel 
$—1 alignment in the ground state, a magnetic structure that 
'^consists of triangles is said to be frustrated In this 
respect the icosidodecahedron belongs to the archetypical 
class of frustrated spin systems made of corner-sharing 
triangles as does the two-dimensional kagome lattice anti- 
Compared to other antiferromag- 
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netically coupled spin systems, such as spin rings for in- 
stance, these structures possess unusual features generated 
by the frustration: (1) many low-lying singlet states below 
the first triplet excitation, (2) an extended plateau of the 
magnetization at one third of the saturation magnetization 



when plotted versus field at low temperatures, and (3) a 
large magnetization jump to saturation again as funtion 
of applied magnetic field. The last feature is intimately 



connected with a huge magnetocaloric effect 15 16 




Fig. 1: The core structure of a Keplerate molecule is an icosido- 
decahedron. The bullets represent the 30 spin sites, the edges 
indicate the 60 exchange interactions. 

For theoretical investigations the extended lattice sys- 
tems such as the kagome excape even a numerical treat- 
ment since the dimension of the related Hilbert space 
grows exponentially like (2s -|- 1)^, where s denotes the 
spin quantum number of individual spins and A'^ the num- 
ber of spins treated in the model. Therefore, the exis- 
tance of finite, i.e. molecular realizations of such frus- 
trated structures provides an opportunity to theoretically 
investigate the quantum energy spectrum and to under- 
stand the related features. A major insight could already 
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be achieved through the study of icosidodceahedra a cou- 
ple of years ago. It turned out that the high-field behavior 
is dominated by special quasi particles, so-called indepen- 
dent magnons 13 17 18 . 

Nevertheless, the full thermodynamics, i.e. physical 
observables as function of both temperature and applied 
field, was so far not available for systems as large as Ke- 
plerates. In this Letter we demonstrate for the first time 
that by means of the finite-temperature Lanczos method 
(FTLM) [19] the thermodynamic properties of the icosido- 
decahedron with 30 spins s — \/2 can indeed be evaluated. 
In an investigation prior to these calculations we demon- 
strated that the FTLM is a very accurate approximation 
scheme that provides quasi-exact results . None of the 
alternative approximations - Density Matrix Renormal- 
ization Group techniques (DMRG) [21 22 or Quantum 
Monte-Carlo (QMC) [23p5) - is able to deliver these re- 
sults for such a molecule. 

Evaluation of magnetic observables. Our nu- 
merical calculations in the Heisenberg model had to be 
performed on a supercomputer. We employed the SGI Al- 
tix 4700 at the German Leibniz Supercomputing Center 
using openMP parallelization with up to 510 cores. The 
complete calculation needed approximately a full week of 
cpu time on 510 cores. The resulting magnetic observables 
are shown in Figures [2j [3j and |4j In case of the magnetic 
susceptibility x we compare to experimental data that are 
published for the highly symmetric molecule {W72V30} 
261. 



interaction J as well as of the spectroscopic splitting factor 
g were adopted. 
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Fig. 2: Susceptibility as function of temperature: Crosses mark 
the experimental data that are published in f26l. The curve is 
the result of our simulation. 

Figure [2] displays the magnetic susceptibility as a func- 
tion of temperature for an applied field of B = 0.5 T. 
The low-temperature part demonstrates that the antifer- 
romagnetically coupled (s — 1/2) icosidodecahedron is 
a gapped spin system, i.e. it possesses an energy gap 
between the singlet ground state and the first excited 
triplet state. QMC calculations cannot resolve the low- 
temperature behavior due to the negative-sign problem for 
frustrated systems. Nevertheless, QMC can accurately de- 
termine the high-temperature behavior as is demonstrated 
in 26 . From this calculation the values of the exchange 
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Fig. 3: Heat capacity: the solid curve predicts the behavior 
of {W72V30}. The dashed curve provides an estimate for the 
phonon contribution. 

The singlet-triplet gap seen in the susceptibility does 
not exclude further singlet states below the triplet. In- 
deed, this is one of the expected frustration hallmarks for 

In Figure [3] we predict that 



the icosidodecahedron 27 



for {W72V30} these singlets contribute dominantly to the 
heat capacity below 10 K. The solid curve, that displays 
the specific heat function at zero field, exhibits pronounced 
Schottky-like peaks which originate dominantly from the 
low-lying singlets. The dashed curve provides a reason- 
able estimate for the additional phonon contribution. In 
an experiment it should be possible to disentangle the two 
below 10 K. The inset shows to full curve up to room 
temperature. 
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Fig. 4: Magnetization for various temperatures: the steps at 
low temperatures result from crossings of Zeeman split energy 
levels. 

Figure [4] provides the theoretical estimates for the mag- 
netization curves at various temperatures. Again the 
singlet-triplet gap is visible, this time as flat M — be- 
havior up to the first magnetization step which happens 
at the level crossing of the lowest {S = 0, M — 0) and 
(S* — 1,M — —1) Zeman levels. In realistic static and 
pulsed magnetic fields the first three steps of the magne- 
tization curve could be observed. 
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Summerizing, by means of the finite-temperature Lanc- 
zos method (FTLM) we are able to accurately evaluate 
all thermal properties of a giant magnetic molecule that 
is structurely related to the kagome lattice. We hope that 
these achivements will promote further investigations in 
the field of quantum magnetism. 

Technical details. — The magnetic properties of the 
Keplerate molecule {W72V30} are described by the follow- 
ing model 



H 



= -2^E 



9^^B 



(1) 



<hj> 



The first term (Heisenberg Hamiltonian) models the 
isotropic exchange interaction between spins centered at 
nearest neighbor sites i and j. J = —57.5 K is the anti- 



ferromagnetic exchange parameter 26 . The second term 
(Zeeman term) represents the interaction with the exter- 
nal magnetic field. The spectroscopic splitting factor is 
taken to be g = 1.95 26 . 



In the finite-temperature Lanczos method 19 the exact 
partition function 



ZiT,B) = EE 



v,r|e-^^k,r; 



(2) 



is approximated by 



ZiT,B) 



Edim('H(r)) 



— 1 71 — 1 



xe-^^""'''|(n(^,r)|z.,r)p 



(3) 



For the evaluation of the right hand side of Eq. ([s]) | i^, T ) 
is taken as the initial vector of a Lanczos iteration. This 
iteration consists of Nj^ Lanczos steps, which span a re- 
spective Krylow space. The Hamiltonian is diagonalized 
in this Krylow space which yields the Nl Lanczos eigen- 
vectors I n(i^, r) ) as well as the associated Lanczos energy 
eigenvalues en'^^ ■ The number of Lanczos steps Nl is a 
parameter of the approximation. Nl ~ 100 yields good 
results [20] . i?r is the number of random vectors that are 
considered in the sum instead of the full basis set. F labels 
the irreducible representations of the employed symmetry 
group. The full Hilbert space is decomposed into mutually 
orthogonal subspaces H(F). An observable would then be 
calculated as 



OiT,B) 



1 



Z{T,B 



E 



dimCH(F)) 
Rr 



EE 



r u=l n=l 

y.{n{v,T)\0\v,T){v,T\n{v,T)) . (4) 



For the present calculations we employed the collective ro- 
tations of the spin about the z-axis as the symmetry. The 
resulting subspaces 'H{M) for the total magnetic quantum 
numbers M — 15, 14, . . . , 0, . . . , — 14, — 15 were included 
exactly for \M\ > 10, since these are small enough. For 
all other subspaces we choose i? = 20 and Nl = 100. 
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